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ABSTRACT 



Elect ionic analog computers have been used ext'snsively for th^ 
solution and display of many dynamics problemSo The majority of the 
problems worked with have been those involving linear differential 
equations with constant coeff icientSo In cases involving non-linear 
differential equations fewer solutions have been developedo This has 
been due mainly to the need of using non-linear elements in the com- 
puter circuits when setting up the solutions. 

One such problem^ that of a mass moving under the action of a 
first power central force, is treated in some detailo The differen- 
tial equation Is derived, the problem is scaled, and the circuitry 
developed. Solutions obtained by the use of the electronic analog 
computer are displayed and compared with solutions obtained by num- 
erical methods and errors and their sources are discussed. Finally 
there is an overall evaluation of the usefulness of analog computers 
to this sort of problem. In an appendix, a second practical dynamics 
problem Is discussed, but a solution was not obtained due to lack of 
time available. 
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lo Tntrodiflct iom 



This thesis presents the electronic analog solution to a non- 
linear dynamics problem which leads to the differential equation 
^ ^ ^ example problem is taken and the steps in reduc- 

ing it to a fom suitable for an electronic analog computer^^ hereafter 
referred to as an analog computer, are showno The results are then 
compared with two solutions obtained by numerical methods^ In Appen- 
dix III an equation of the form; X ^ A X f 8 A f h ^ O 
discussed and the problems encountered in trying to obtain an analog 
computer solution are delineatedo 

The writer wishes to express his appreciation for the assistance 
given him by Professor John E* Brock and to the Professors, particular- 
ly Professor 0<> H*. Polk, and the technicians of the Electrical Engineer- 
ing Department o The numerical solutions in Appendix I were contributed 
by Professor Brocko 
2o Background, 

Solutions for many dynamics problems have been established using 
the analog computer and references can be found in the technical litera- 
ture, One such reference for a non-linear problem. Analog Computer Solu- 

2 

tion of a Non-Linear Differential Equation, by H, G, Markey, ( 2) , was 
found but was only applicable in a general way to this invest Igat ion. 

It was considered that if a means could be found to display some 
of the classical problems encountered in early college dynamics on the 
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2 d/dt Is denoted by a dot placed above the variable operated on. 
Numbers In parentheses refer to references in Bibliography, 
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ajinaiofe. computer the following advantages would be obtained; 

(a) the general usefulness of the analog computer could be 
made more readily apparent; 

^b) in dealing with these problems attention could be focused 
on the dynamic principles leading to the governing differential equations 
and upon the mechanical significance of the results and not upon the mathe- 
matical difficulties In obtaining an analytical solution; 

fc> in the case of those problems where analytical solutions 
have been obtained for certain particular parameter values and which 
therefore seem to be separated into many different cases the dynamical 
significance of which is not apparent, the general problem could be 
dealt with directly; 

fd> it would be possible to include normal dynamical in- 
fluences fsuch as energy loss due to pivot friction or windage'^ ^17ich- 
out so complicating the mathematics of the solution that the modified 
problem appears to be entirely different from the Idealized problemo 
In addition to the above it was desired to obtain these results 
using only the analog computers and their associated equipment normal- 
ly available In an analog computer laboratoryo 
3o The problem and general method cf solution* 

The problem considered was that of determining the subsequent 
motion of a body weighing 1930 pounds, attached to a spring having a 
free length of five Inches and a scale of ten pounds per loch, when 
released with the following initial conditionSo At the Initial instant 
the radius is four Inches and its rate of change is zero; the polar 
angle, 0 , is zero and Its rate of change is three radians per second* 
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TH- - 7 'tlter end of the spring Is attached to a flKed point and the body 
IS ^v6rmitted to slide v^ithoiat friction upon a horizontal planeo 

We will discuss the sequence of steps necessary for the solution 
oi a problem of this type^ and then we will proceed with the solution,) 

One might of course proceed with a foil scale experimental program as 
a method of solution^, but eliminating this possibility we would^ 

Be Using the p>rincipals of Mechanics arrive at one or more 
differential equations describing the motlono 
bo Solve these equations^ incorporating the starting condi- 
tionso This solution may be analytic^ numerical^ experi- 
mental (dealing withj possibly ^ scaled down mechanical 
variables)^ or by means of an analogy In which one deals 
experimentally with variables of another type (such as 
electrical) which satisfy similar differential equations^ 

Co Interpret the mathemat ical^ experimental ^ or analog results 
in the proper mechanical light so as to arrive at a mean- 
ingful solution to the original problem^ 

In this thesiSg we are investigating the practicability of pro- 
ceeding immediately from the first step to a solution by use of stand- 
ard analog computer equipmento We do not have an analytical solution 
to the problem stated^, although it is likely that one might be obtained 
in terms of elliptic functions and integrals, Howeverj in an appendix 
we will develop two different numerical solutions to the problem with 
which we can compare our analog solutions , 
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The differential equation of motiono 



Figo 1 shows the body in a general position^ The solid arrow re- 
presents the spring force = 10(r-5), where r represents the radial 

distance from the fixed point Oo The dotted arrows represent D'Alembert 

forces necessary for dynamic equilibriumo We see that 

-r yvi ::: O 

YT) = O 

Now by kinematics, ^ and , From the 

' /7 ciit* 

second equation we see that ri '^0 ~ C = Cc'/iiZi, « This can also be seen 
from the fact that the angular momentum of the system about 0, namely 
, is invariable. From the first equation, we have 
lO j-jO , and from this we get '2h—l0 t/7 ~hO ^0 . 

S % ' lr 

<? V ^ ‘ L. 

Substituting O’ - we finally get —^■'7 i~ /O 

nx 

In our case, evaluating C at the initial instant we have C = 48. Thus 

we have as our set of differential equations: 

A - -A 3,c y ^/o 

Now it is possible to perform some mathematical manipulations which 
simplify this system. In particular, it is easily possible to obtain a 
first integral of the first equation of the system. However, we regard 
it as contrary to the spirit of this thesis, the purpose of which deals 
with the ready feasibility of making an analog computer solution of this 
system, to perform any such manipulations, and it Is 
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Fig. 1 Force Diagram for Central Force Problem 
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this system with which we shall be directly concerned when we attempt 
the analog solution* The numerical solutions for this problem will be 
found in Appendix I. 

5^ Discussion of equipment* 

Before taking up the solution of the problem, a description of 
the equipment used in the solution of this problem will be presented* 

It is assumed that the reader is already acquainted with the basic theory 
of the analog computer and with the usual circuitry used, such as summers, 
integrators, etc* The Handbook of Automation, Compution and Control, 

Vol* 2, E* M* Grabbe, (1), is a good reference for this information as 
well as for additional information on the equipment discussed below* 

A* Donner Analog Computer, Model 3000* 

This analog computer, which can be used for the quantita- 
tive solution of linear (and certain classes of non- linear) differential 
equations and transfer functions, contains ten DC operational amplifiers, 
any one of which can be used for addition, subtraction, multiplication or 
division by a constant, sign changing, or integration. Problems express- 
ed as differential equations are entered in terms of electrical components 
on a detachable problem board. Stability and accuracy are satisfactory 
for problem solution times up to 100 seconds or more which permits ac- 
curate recording with conventional pen recorders. (5) 

B* Donner Function Multiplier, Model 3730* 

This function multiplier consists of two multiplier 
channels and a regulated power supply* Each multiplier channel produces 
an output voltage which is accurately proportional to the instantaneous 
produce of two independent input voltages, where each input is either 
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ceostant or varying with timeo Either input may be positive or negative^ 

5 ;-:^ that four quadrant multiplication is providedo The range of output 
and input voltages is plus and minus 100 volts; this being necessary to 
stay within the linear range of the operational amplifiers of the com- 
puteto To maintain the output voltage at 100 volts or less the Function 
Multiplier is designed to give an output voltage which is <,01 the product 
of the input voltagaSo ( 6 ) 

For the solution of the problem of this thesis two of these multi- 
pliers were usedo They gave accurate results when used for straight 
multiplication although they do tend to drift over a period of time and 
have to be readjusted; this is a minor operation^ however <, l-Jhen used in a 
division circuit^^ which is discussed in a later sectionj> the results obtain 
ed were not as accurate^ however o It is believed, however^ that this was a 
fault of the circuit and not of the ^unction multiplier bacausa^ as mention 
ed abovej> the function multipliers gave quite accurate outputs when used 
for multiplication alonCo 

Co Donner Function Generator, Model 3750« 

This variable base function generator is designed for use in 
conjunction with two external operational amplifiers to produce an output 
voltage which can be adjusted to approximate a desired single valued 
function of the input voltageo One operational amplifier is required for 
operation of the function generator and the other is used to accept the 
output sigml at its summing junction<. This amplifier may also be used 
for additional summing, inverting, integrating or other operationSo The 
function generator operates on the principal that the function can be 
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approximated by a series of straight line segmentSj^ each segment being 
limited to a slope between plus and minus two volts per volto The Input 
and output voltages may vary between plus and minus 100 volts« (4) 



For the solution of the thesis problem it was desired to use 

3 

this function generator to generate the function 2304/r but it was found 
that the slope of curve for this function exceeded the capability of the 
equipment. This is duscussed further in Sec. 9, If it had been possible 
to use this function generator the two function multipliers would not have 
been required, 

6, Computer equation and scaling. 

To reduce our problem to a form suitable for the computer it is 
necessary to apply scaling factors. This was done using the methods out- 
lined in Basic Theory of the Electronic Analog Computer, by R, C, H, l^eel- 
er^ (9), A brief summary of this process is presented here. 

The differential equation to be solved is first arranged so as to 
solve for the derivative of the highest order. In our case ^ 

-BX + C . The equation is then scaled so that maximum value of each para- 
meter is represented by a voltage, close to but not exceeding 100 volts. 

To do this scaling factors are assigned as shown by the following ex- 
ample: 

X ^ X 



Here X is the computer voltage representing the variable x, and 
is its scaling factor. After being calculated the scaling factor is 
usually rounded off to facilitate computations. After suitable scaling 
factors are found, the equation is put into the form: 
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To develop the applicable circuits for the problem solution it is 
necessary to determine the values of resistance and capacity needed 
for each component of the circuit. Using the procedures in iJheeler®s 
bookj f9), pages 2-10 we find, for example, that an operational amplifier 
when used as a summer has an output voltage ) 

or in our case r: — /I- ^ • If we now let Cc^ -- 

where a^ is a coefficient potentiometer setting and and are re> 



sistances, we can establish the relationship 






It should be noted here that an R with a subscript, R^* refers to a resis- 
tor and without the subscript refers to the voltage representing the 
variable r, the radius of the problem. Now the above relationship can 



be solved for a . For integrators the relationship is - 






where refers to a capacitor. 

7« Analog computer circuits. 

In Appendix II the calculations for scaling the differential equa- 
tions of our problem are presented. After scaling we have the follow- 



ing equations: 



= B- i-2a 
d- 3o,aooy^<_ 



. S y L'jOoO 

K5“ 



Before solving this problem on the analog computer two main decisions 

3 2 

have to be made: first how to calculate R and R , and second how 

3 ^ 2 

to develop the terms Z = 576,000/R and & = 30,000/R . It was 

hoped at first that the terms for Z and 3 could be developed using func- 
tion generators but as mentioned previously this proved unsatisfactory. 
Thus it was expedient to use the division circuit shown in Fig. 2. 
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Division Circuit (6) 
Fig. 2 



With this circuit Z = 



lOORj X 



The factor 100 results from the 



R Y 

2 3 

output of the multiplier being .01 XI. If we now let Y = R and X 
= constant, using the above relationship we should be able to develop 



Z = 576,000/R . 



We know from the parameters of the problem that when r is 4^ 

R should be 20 volts. If we then put this value through two function 

2 3 

multipliers we come up with K R . As this value is small, ,8 volt, we 
multiply it by a factor of two using an operational amplifier and then 
put it into the function multiplier of the division circuit. Also using 
this value of the voltage for R we can calculate the value Z should 
have, in this case 72 volts. With these values we can now solve for a 
value of X so that with an input of 1.6 volts for Y and the calculated 

value of X , Z will be 72 volts. Solving for X: 

X ^ ^ ^ 

/ ~> fKf ■ i\, /oO fi, 

Now letting R ^ equal lOM and R equal IM, wa find that an X of 

11 o 5 should be used. (It was found that resistances of lOM and IK i^orked 

better than resistors of O.IM and IM) . This same procedure was applied 
to 6^ and the corresnonding voltage, X^ was found to be 30 volts. 
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It should be noted here that another method for determining B 
presents itself, that of multiplying Z by R/19o2. By doing this the 
second division circuit could be eliminated and only another multiplica- 
tion, with its more accurate results, requiredo This method was tried 
and it was found that for some unexplained reason 0 passed through 
zero and became slightly negativee As a result of this ^ oscillated 
instead of increasing smoothly from zero to a maximum value o For this 
reason the division circuit for developing & was used. 

After the above determinations were made, the circuit of Fig« 3 was 
assembled and computations made* In assembling the circuit the values of 
the a*s calculated in Appendix II were adjusted for the actual values of 
resistances and capacitors used, e«go, 1.005 M actual resistance vs. 
nominal value of IM. Figs. 4 and 5 are photographs of the setup used 
and shows the relative simplicity of the final setup for solution of the 
problem. 

8. Results. 

After assembling the circuit of Fig. 3, it was found that to obtain 

the desired values of voltage for Z and ^ the values of the input vol- 
tages calculated for X and had to be adjusted. For X, a value of 

20 volts, and for X^ a value of 33 volts was required. Once these ad- 

justments were made, the computing runs were made and the results are 
shovm on the Brush Recordings of Fig. 6 and Fig. 7. These recordings were 
all made using a paper speed of 5/mm/ sec and with varying voltage scales as 
shown on each trace. These results are also summarized in the table of Fig. 
8. From these results curves were plotted and then compared with the re- 
sults obtained by the numerical solutions^ as shown in Figs. 9, 10, & 11, 
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In analysing the results each term will be considered separate- 

4 *> 

ly. Considering r first it is seen that the maximam value of 38 ob- 
tained agrees with the maximum value of the numerical solution but that 
the minimum value of - 12.5 is lower than the - 15.26 of the numerical 
solution. This latter discrepancy is attributed to the actual values 
obtained for Z and will be discussed later. 



11 




f 












A or. 




12 



Fig* 3 Circuit Diagram 
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Fig 11, Radius vs Angle of Rotation 
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Fo'r r the value obtalried was 9o7 and the ml) .unuiB va.ue'fe 



w-^i-re zero when r was a maximum and a mlnimumo This agrees veil with 
the numerical solution where the maximum value was 9«6c Considering r 
we see from Figo 9 that the analog values are slightly higher at all values 
than the r®s of the numerical solutioBo This error is not considered 
excessive o 

The largest discrepancies appear when we consider ©o As can be 
seen In FigSo 85 IO 5 and 11 the analog value reached its maximum for 
the first apse (point of greatest distance from the center of atf^raction) 
rapidly and then remained constant for a period of tim€o Here as with 
r the discrepancies are considered to be caused by the values obtained 
for Go 

Considering the problem overall^ the more significant results ob- 
tained appear to agree rather well with the values obtained by the 
numerical solotionso The major discrepancies appear when the part of 
the circuit handling the division is consideredo As can be seen from 

c 

FigSo 6 and 12 for and Figo 7 for the outputs of these division 
circuits change rapidly to a small negative voltage and than remain re- 
latively constant for a period of timeo We can also see from Figo 10 
that the division circuit does not do what theoretically it shouldo Thus 
for either parameter the minimum voltage desired^ when r is a maximurf^it 
is never obtainedo With Z^ this term is small when compared with the 
others in summing for r and the effect is not pronouneedo With © 
however this defect has a more pronounced effect and © is not developed 
in the smooth curve desiredo 
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9 Discussion of discrepancies o 

The discrepancies found in the above problem wen- a£t£ib^p'-“ 

ed to the division circuits usedo No satisfactory answer could be found 
as to why the desired divisions could not be obtained « It is known 
that a circuit such as this develops a certain amount of noise o That is^, 
the function multiplier has a certain amount of noise inherent in its 
output and that if this is put through an operational amplifier this 
noise is amplifiedo The Handbook of Automation^, Computation and Control^ 
Volo by Grabbe fl) discusses this briefly and mentiors that a small 
capacitor placed in parallel with the multiplier will help to alleviate 
this problemo This was tried but did not give satisfactory result So 

As mentioned previously^ if a function generator could have been 
usedj the circuitry could have been simpllfieds, no function multip- 

liers would have been required o With the Donner function generator the 

3 

slope of the function 576s, 000/ R j, for low values of Rj exceed the maxi- 
mum of two volts per volt permitted by the device® One other type of 
function generator was tried® This was an Autograff Tf plotter converted 
to a function generator by replacing the recording pen with a' pick-up 
coil and plotting the desired function with a conducting Ink® However^ 
with this arrangement the desired range of voltages could not be obtain- 
ed® 

Still another type of function generator that might have proved 
sat Isfactorys, if it had been available^ is the photo-former type® This 
type of function generator operates as follows® The basic piece o£ equip- 
ment is a cathode-ray tube® An input voltage is applied between the 
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i. IfSctiom plates c£ a cathode ray tube thr:>ug’ suitabi^^ 

The voltage between the vertical deflecti- platen la 

as the output voltageo This voltage is made to vary .is a fanctio! 
of the input voltage by a feed-back arrangement which fotces the electron 
beam to follow the boundary of an opaque mask placed over the lower poi- 
tion of the cathode-ray screeBo Thus as the spot on the cathode-ray tube 
screen emerges from behind the mask a photocell in front of the tube ap- 
plies an error voltage across the input terminals of the vertical dafiec- 
tion d-c amplifier 5 so phased that the beam is forced downward toward the 
face of the masko Therefore if the mask is shaped to represent the function 
being generated the spot will follow this curve and deliver an outpvit 
voltage proportional to the input voitageo This type of function genera- 
tor ioS said to be very accurate in developing many functionSo (3) 

10« ConclusionSo 

Considering the results obtained from this problem (keeping in mind 
that indeed it is but a single problem) ^ it was found that a typical** non- 
linear dynamics problem can be set up on an analog computero However this 
type of set-up is not done rapidly or easilyo Considerable thought has 
to be given as to what type of equipment shall be used and what kind of 
circuits are nacessaryo Because they require the use ©f various types 
of non-linear computer accessories the circuits become very sensitive and 
results accurate to the degree normally expected from the analog computer 
may not be obtainedo Care has to be taken in selecting scaling factors^, 
where powers and roots are involved^ to avoid over- loading the operation- 
al amplifierSo It was found^ however g that the function multipliers used 
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c R gave quite, accurate results^) sve^ni at ly^- •■-It- 
they ‘were kept balancedo 

la setting up a problem of this type it will usually be found 

3 

that tfere will be one key term to be developed^, such as the A/R of 
this problemo Once a way is found to develop or represent this term 
the remaining computer setup is routine and with patience and luck a 
solution can be obtainedo 
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APPENDIX I 



Nv.metlcal lutjLons for Central Force Problem 

The atatement of the problem is given in SeCo 3 on page 2« Re- 
presenting this Inlonnat ion in mathematical terms, we have F , 

^ ® ^ spring 

2 

10(r-5) IbSo and m • 1930/386 s 5 lbs sec /in^j, and initially (at 

a r 

time = we have r = 4 incheSj r = Oj, 0 ^ 0^ and Q ^ 3 

radians/seCo Since energy is conserved, we have 
g _ j ^ ^ ^ r (^7 -s ') - C rJ 

where E, T, and V are expressed as energy per unit mass in units 
of Ir^/sec^o Here we have used \/ - i - •' 

Using the initial conditions to evaluate E, we have 

L ‘AA M'.i f ■' - •> 

. "V 'X / 

... '•/ ] 

i 

Ac apses ^ Jo.- Apsidal radii are 

j. \'t 

given by f j - 7/i Substituting 

~ /' /f and rationalizing we get: 

*/• - ';.3 

; - y, •) ^ X!^rff -.t. ; ^ 

- / « •'/ ^ ' W ^ , 

K'tiiowing one: root i,s 4 we obtain (/7 — J ( h • . : ■* /> .' 5 ;,. .r 

This can have only one positive roota Synthetic division indicates a 
root of appjroxiroately 13o2 and using Newton® s method: 



2: 



1 

1 



r 



t 



- i^n -'yi. 

o^ , - / -"J ■“ 3- O 3 ) . ;,-i ._ C~'^') 

'C^'C/3) 

cx._j_ ^ 13 ./^' - 3 ^ 63 ./>) , y3.o~'_ //3~y^ _ /3,/vy3'" 

^srr 9 > 

cx^ -. 13 ./V^/r_^((3./yyj-J /)./v7r- ■ O 

/vyj)':) ~ 

- ;3 ; /y yo 3 

To find the apsidal angle and r and 0 as functions of time^, we 
resort to a numerical procedure since the integral involved is not 
elementary o Returning to fundamentals we have: 

"T-hV J t- 0y~s-)~^~'?3 

K - - 9^ 

/.• P - 2xA£^ ~h/0 

/I ^ 

We also know = 4 and = 13«144. Now using an iterated 

integration, a curve of r = r(t) is assumed such that r 0 at 
the end points (apses><, The apsidal time is divided into n equal 
intervals /r\ ; c being as yet unknown« We will use n = 6, although 

a larger n will give a more accurate result 0 

Assumed values of v® are selected for each epochs Values of r 
are calculated and integrated with the condition at ^ 

This should yield/7'^? at 6-^9 but there is an error ^ 0 We remove 

this error by using a correction curve which is essentially A ^ z 

expressed however in appropriate form for and obtained by numerical 

«» 

integratioHc This arises from assuming that the error in r is due 

Cf v/ 

to an error in r which must be essentially parabolic in naturej vanish- 
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trig at the end points since the apsidal distances are known The rest 
of the calculations are self explanatory and lead to the curved path 
shown in Figo 9<, 
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Second Numerical Method 



From the expression for 1 ^ 

Upon substituting we get as a function of t 

and thus can construct a curve of /7 as a function of (We take 

the positive branch of the square root so as to deal with the period 
during which r is increasing from 4 to 13.145 inches.). Also we have 
=. i~ / 0 so that we can construct a curve of a/ 

iC 

as a function of r, and this relation shows that // — C when r is 

^ 6 C 

approximately equal to 7.62. Having curves of both // and ; as 

C 

functions of r, we can construct a curve of h as a function of A 

The differential of time may be written in either of two w^ys 

i)^ i) A 

9 and this permits us to write 

A /T 




so as to avoid infinite values for the integrands. These calcula- 
tions can be carried out by numerical methods as shown on the follow- 
ing pages and ^ is found to be 1.4863 seconds j, which agrees with the 
value of "h:. calculated in the first numerical solution. 
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APPENDIX II 



Scaling Equations for Central Force Problem 
The basic equations for this problem are: 

(9 -- 

The initial conditions at time t = 0 are: 

/T-Vj ^ ^ O ^ <3 --0 ^ 

and we also know from the numerical solution that the approximate 

maximum values of each of the parameters are: 

A7 - /3 A 

Jl - 3 9-3 

Also from the numerical solution we know that the time from apse to apse 

is about 1.46 seconds. For our solution we will select (1.5 

revolutions) as we are only interested in the initial aspect. Knowing 

the maximum values and using the relationship x 3 as explained 

in SeCo 6, we can now solve for the scaling factors. 

/, Le.9 - I X V/ 

/, L’e.ir Z ,/Z. 

/, A" Ay,^ ' 

' ^ l_c ^ oC^ - 2. ;)£j i/ 

/ , L.C.L ^ . 0 </ 

problem time down so we assume 
computer of 7.3 seconds^ apse to 





- n.A 

/ o o 








. 09 - 




/o o 






- _3jr 


^ , 3 r 


/ OC 






f OO 


- , 09 / 1 . 




- ^ 
/c^o 


.0 3 



For time scaling we wish to slow the 
- 5"^^ s giving a period for the 
apse. 
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Now using the above scaling factors we can proceed to scale 



the basic equations for the computer: 



R 


_ '220^ 


<loC^ R. 




-h 


/ 0 












i. 

n 


~ a3cV 








-B. 














_ 2-0.2. 






/6 










f- 






■3 ^')(^jOOO 






■h 


'3.0 



S -- 

- - 

yO‘/C^f-'(V- 

- 'id doo 



The equations are now in the form ~ '--(uj. Sr CjU^ ^ . 

where Z = 576,000/R , and = SO^GOO/R , Because Z and ^ are 
developed by the division circuits and 20 is a constant voltage their a^ 
values (coefficient potentiometer settings) are each lo Thus their 
corresponding resistors are all IMo To detemine the value of we 
use the relationship • Equating to 

and letting R 2 and equal IM, we find to be .8. 
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To obtain R and R we must integrate ^ f 
vJe therefore scale these equations as follows; 

■X ,; k - 

-t- 

r A i^c 

AcJL, J 

n r. Ji .^cj^^a<K 

^/, - ~S ^ 

n^o 



oCj 

A - ^ 

^4- 

- _ / b' 






-1^ \ R 



-ii \aM^---,Ii-\p-.' 

. I ' 






Now using the relationship - 



c;t 






we solve for a. and a_ 

H ! 



as shown below* The resulting values are *833 and *48 respectively 
using the resistors and capacitors shown* 

To complete our scaling we must now determine tl e a^ value 
for ^ is found by integrating ^ * Scaling this 



equation 



we finds 







I 

I 

1 

I 



3 



OJ ^ S' 3 3 

^ S' - ~ 



OUcf - 

«, Cf ■ 





CT ^ ^ 














a., -- ^'>1 





,'. a, ^ ,30 



The problem is now ready for the computer« The table on the 
following page summarizes the values for all resistors and capacitors'* 
used for this problemo 
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Table of C Ire jit Elements 



Amplifier 
(See Eig, 3 ) 


Circuit Element 
^nd Value 


! 

Remarks 


1 


R7I = IM 
R72 “ 10^^ 


Forms *=Z 


2, '• 


R]_ r IM 
Rg : iM 

R^ - IM 

R^2 " 

= 1 0 0 

02 “ 0 0 8 

03 s 1.0 


• * 

S urns + R 


3 


R|^ • - IM 
C ^ 2 ~ 1 jLi-i* 

9^1 - 0.833 


- J" +R dtp = -R 


h 


R^ = 2 M 
Cfi; " 2;jf 
0^ “ Ooij.8 


- J -R dtp = +R 


5 


Rp = 5 M 
= ljUf 

0 Q — 0.33 


dtp= e 


7 


RqI “ IM 
R32 - lOM 


Forms — ^ 


9 


R^ = O.IM 

R - 0 . 2 M 

f 9 

05 - 1.0 


2 r 2 r 3 



Fig. A21 
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APPENDIX III 



The Spinning Top 

During this invest igation^ some attention was also directed to 
the problem of the spinning top; that is a symmetrical top with one 
pc int fixed under the action of gravity. Using the notation of Synge 
and Griffith, (7), the equation ^ (/ -j^) 

was reduced to the form of 

A sample problem was then selected from Applied Mechanics by No C, 
Riggs^(8), which consisted of a gyroscope being released with & equal to 
60 degrees and the subsequent motion being an oscillation between 60 and 
82 degrees while precessing at a variable rateo 

With this problem^ as with the central force problem^ the powers 

of X could be obtained using function multipliers. The key to the 

^ 2 

solution however consists of developing the square root of X <> This 
should be possible using either a function generator or a division cir- 
cuit similar to that of the central force problem. Using the function 
generator was tried first, A curve of Y = 10 X 1/2^ X and Y 

being arbitrary voltages from zero to lOOj was calculated and set into 
the function generator. It was then found on scaling the problem that 

to prevent the term CX from exceeding 100 volts it was necessary to 

# 2 

assign to X a maximum voltage of approximately nine volts. Then 
when the circuit was assembled it was found that this voltage was too 
small for the function generator to operate. 

When the use of the function generator proved to be unsatisfactory 
the use of a division circuit was attempted. To use this it was assumed 
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that the following relationship was valid. We know that for this 

circuit Z = J_0C Now letting Y equal Z and then solving 

^ 2 ^ 2 

for X we find that X = R-> Z , or Z equals 10 X * ^ with 

lOOR^ 

s R^ = IM. This arrangement was then put into the circuit and 
computations attempted. It was found that the operational amplifier 
generating X went from its initial value through zero and then over- 
loaded, i.e., the voltage X (representing x = Cost ^ ), did not 
vary as it should since the circuit was unstable. 

After the above attempts failed it was decided that this particular 
problem could not be scaled satisfactorily for an acceptable analog 
computer solution within the time available. 
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